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Abstract 

In this note we propose a transformation which decouples stationary Mean Field 
Games systems with superlinear Hamiltonians of the form |p|^ , r' > 1, and turns 
the Hamilton-Jacobi-Bellman equation into a quasi-linear equation involving the 
r-Laplace operator. Such a transformation requires an assumption on solutions of 
the system, which is satisfied for example in space dimension one or if solutions are 
radial. 
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1 Introduction 


Mean Field Games (briefly MFG) is a branch of Dynamic Games which has 
been proposed independently by Lasry, Lions 0. i. 0. mi and Gaines, 
Hnang, Malhame [6], and aims at modeling and analyzing decision processes 
involving a very large nnmber of indistingnishable rational agents. In MFG, 
every agent belonging to a popnlation of inhnite individnals has the goal of 
minimizing some cost which depends on his own state and on the average 
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distribution of the other players. Suppose that the state of a typical player is 
driven by the stochastic differential equation 

dXg = —agds + y/^dBg G Q, 

where as is the control, Bg is a Brownian motion, z/ > 0 and the domain C 
M'^, d > 1, is the so-called state space. Suppose also that the cost functional 
has the long-time-average form 

jr(Xo,a) = liminf [ E[L{as) + fiXs,ms)]ds, 

T^oo i Jo 

where the (convex) Lagrangian function L{a) is associated to the cost paid by 
the player to change his own state, and the term involving /, that we assume 
to be a X [0, oo)) function, is the cost paid for being at state x ^ fl, 

and it depends on the empirical density of the other players. Then, under 
the assumption that players are indistinguishable, it has been shown that 
equilibria of the game (in the sense of Nash) are captured by the following 
system of non-linear elliptic equations 

/ 

(HJB) —vXu{x)^B{Du{x))^\ = f[x^m{x)) in 
< (K) —i'Xm{x)—dw{DH{Du{x))m{x)) = Q in (1) 

J^m{x)dx = 1, m > 0. 

Here, H denotes the Legendre transform of L. The two unknowns u, A in the 
Hamilton-Jacobi-Bellman equation (IT])-(HJB) provide respectively the optimal 
control of a typical player, given in feedback form by a* : x i—)■ —DH{Du{x)), 
and the average cost J'{Xo,a*). On the other hand, the solution m of the 
Kolmogorov equation ([T])-(K) is the stationary distribution of players imple¬ 
menting the optimal strategy, that is the long time behavior of the whole 
population playing in an optimal way. The two equations are coupled via the 
cost function /. Note that a set of boundary conditions is usually associated to 
o. for example u, m can be assumed to be periodic if = (0,1)'^ or Neumann 
conditions are imposed if is bounded and Xg is subject to reflection at dfl. 
In some models is the whole M'^. 

A relevant class of MFG models assumes that the Lagrangian function L has 
the form L{a) = /q > 0,r > 1. Consequently, the Hamiltonian function 

H becomes 

H(p) = ^ipr', fto = ly > 0, / = -C > 1. (2) 

In the particular situation where L and H are quadratic (namely r = r' = 2), 
it has been pointed out (see 0. n) that the so-called Hopf-Cole transfor¬ 
mation decouples ([I]), and reduces it to a single elliptic semilinear equation 
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of generalized Hartree type. Precisely, let (p := ce 2 “ , c > 0, then (IT])-(HJB) 
reads (setting for simplicity ho = 1 ) 

— 2u'^Aip + {f{x,m) — X)(p = 0 in n (3) 

for all c > 0. Moreover, if we set = (^/q , an easy computation shows 

that is also a solution of (IT])-(K), so if uniqueness of solutions for such 
equation holds (that is true, for example, if suitable boundary conditions are 
imposed), then m = (^^, and therefore cp becomes the only unknown in ([3]). 
This transformation can be exploited to study quadratic MFG systems, both 
from the theoretical and numerical point of view. This strategy is adopted, for 
example, in the works n, 0 , a, 0 . 

The aim of this note is to show that if r' 7 ^ 2 there exists a similar change of 
variables, involving a suitable power of m, that in some cases decouples ([ 1 ]) 
and turns the Hamilton-Jacobi-Bellman equation (IT])-(HJB) into a quasi-linear 
equation of the form 


I-/lArip + (f(x,cp’') - A)cp’' ^ = 0 in Q, 

\f^cp^dx = l, cp>0, = \ 

where Ar(p = (\\Y{\DipY~‘^Dip) is the standard r-Laplace operator (r is the 
conjugate exponent of r')- Such a transformation is in particular possible if 
the vector held vDm+DH{Du)m, which is divergence-free because of ([T])-(K), 
is identically zero on hi. 

Let us briehy recall in which sense {u,m, A) solves ([ 1 ]). 

Definition 1.1 We say that a triple {u,m,X) G (7^(11) x iT[Qj,^(r2) x M zs a 
(local) solution of ([I]) if u, X solve pointwise (II])-(HJB) and m solves ([T])-(K) 
in the weak sense, namely 

V ( Dm-Di+ f mDH(Du) -0^ = 0 G (5) 

Jn Jn 

We say that a couple {ip, X) G (iT[Q^(r2) fl L“^(n)) x R is a solution of (jl]) if 

p / iDipf-^Dip • + / {fix, ip^ - X)ip^-^^ = 0 Ve G C^{n). 

Jn Jn 

Then, the transformation can be stated as follows. 

Theorem 1.2 Suppose that H satisfies fl2l) . 

a) Let {u,m,X) he a solution of ([I]). If the following equality holds, 

uDm + hom\DuY ~‘^Du = 0 a.e.inVt, (6) 
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then 


( 7 ) 


if := in Vt 

is a solution of (jl]). 
b) Let {(p,X) be a solution of (jl]), and suppose that there exists u G 
such that 

hQ(p\Duf ~‘^Du + urD(p = 0 inVt. (8) 

Then, {u,m, A) is a solution of ([T]), where m := in fl. 

The proposed transformation reveals a connection between some (stationary) 
MFG systems with non-quadratic Hamiltonians of the form ([2]) and r-Laplace 
equations, which have been widely studied in the literature and appear in 
many other areas of interest. Apart from existence and uniqueness issues, 
this link might shed some light on MFG problems in general, which can be 
translated into problems involving the r-Laplacian (e.g. qualitative properties 
of solutions, MFG in unbounded domains, vanishing viscosity limit z/ —)■ 0). 

In the quadratic case r' = 2, we have mentioned that if solutions of ffTl)-(K) 

Hqu / hQU \ ^ _ 

are unique, then m = e“ 1 .In this case condition ([ 6 ]) is easily 

verihed, and the assertion of Theorem 11.21 is that (p = solves (jlj) with 
r = 2, which is precisely ([3]). In this sense our transformation can be seen as 
a generalization of the standard Hopf-Gole. 

Note that the change of variables m = ip^ is local, in particular it is independent 
of boundary conditions that might be added to ([T]). However, in order to verify 
dfi]), dH) and therefore to apply Theorem ll.21 it is necessary to specify additional 
information on the problem. Space dimension d = 1 and Neumann conditions 
at the boundary, or u, m, (p enjoying radial symmetry are two possible scenarios 
where (| 6 ]), (jS]) hold. 

Corollary 1.3 Suppose that H satisfies ([2]) and G = {x G : |x| < R} 
for some R G (0, cx)]. Then, {u,m, A) is a radial solution of ([T]) if and only if 
{ip, A) is a radial solution of (jl]), where p = 

Corollary 1.4 Suppose that H satisfies ([2]) and G = (a, 6 ) for some —oo < 
a < b < oo. Then, {u,m, A) G x x M is a solution of ([I]) satis¬ 

fying the Neumann boundary conditionl^ u'{a) = u'{h) = m'(a) = m'(b) = 0 
if and only if (p, X) G fl is a solution of (jl]) satisfying the 

Neumann boundary conditions p'{a) = p'{h) =t), where p = . 

Acknowledgements. The author wishes to express his gratitude to the ano¬ 
nymous referee for his very careful reading of the manuscript and his valuable 
advices. 

2. For equations ([I|)-(K) and ([3|), Neumann boundary conditions are intended in 
the weak sense, namely the space of test functions ^ is set to be 
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2 Proofs 


Remark 1 We point out that m G (f2) for all g > 1 by standard regularity 
results on weak solutions of Kolmogorov equations. Moreover, the Harnack 
inequality guarantees that m > 0 on However, even if m G we do 

not expect in general the same regularity for m, since DH{Du) might lack of 
the desired smoothness if 1 < r' < 2. If r' > 2, it is possible to conclude that 
m is twice differentiable on H and it solves ([T])-(K) in the classical sense. 

On the other hand, it is known that a solution g? of (jl]) enjoys local 
regularity (see, for example, 0. CDl). 

Proof of Theorem II.21 a). Note that m G for all g > 1 (see Remark 

H]), and therefore has the same regularity. Moreover, m > 0 on O, hence 
(p > 0 as well. Equalities (jb]) and ([7]) imply that 

ur—— = —holDul^ ~^Du a.e. in O. (9) 

We multiply the Hamilton-Jacobi-Bellman equation ([T])-(HJB) by where 

^ G 0^(0) is a generic test function. Integrating by parts, 

i, f Du- Dfe'-*) + / f^|Z)«(x)rV-‘ + / (A - = 0. (10) 

Jn Jn r Jn 

We note that by ([9]) 

z/ / ^Du ■ D{ip^-^) = u{r - 1) [ ^Du ■ Dip f ^\Dufp^-\ 

Jn Jn r Jn 

and (r')“^ = (r — 1)^“^, so (fTOj) becomes 

V j Du-Di + / (A - = 0. (11) 

Jn Jn 

Again using ([9]), we obtain 


/ \DpY-‘^Dp . 

Jn yur J Jn 

= !(-) /(A-W, (12) 

v \i>r j Jn 

since (r' — l)(r — 2) + r' — 2 = 0, and by virtue of (ITTD . Equality flT^ . which 
holds for all test functions is precisely the weak formulation of (jl]), hence we 
are done. We hnally observe that p enjoys local regularity (see Remark 
H]), which is inherited by m. 
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b). It is easily verified, in view of ([H]), that m solves (IT])-(K). Since ip = m}!'^ 
is positive (see Remark [T]), 

= —hQ\DuY'~^Du in Vt. 

By carrying out backwards computations flTOll - ffT^ of part a), it follows that 
M, A is a weak solution of ([I])-(HJB). Standard regularity results for the Poisson 
equation guarantee that u G C^(r2) and ([I])-(HJB) is satished in the classical 
sense. □ 

Proof of Corollary 11.31 In the following, p := |a;| and = ep{x) = x/p will 
denote the standard unit vector in the radial direction. 

Let {u,m, A) be a radial solution of ([T]). We write u{x) = u{p), m[x) = rn{p), 
Du{x) = u\p)ep, Dm{x) = m'{p)ep. The Kolmogorov equation ([I])-(K) reads 

/ {uDm + ho\DuY'-^Dum) -0^ = 0 (13) 

Jo. 

for all test functions ^ G (^“(fl), and 

i>Dm{x) + ho\Du{x)Y ~‘^Du{x) m{x) = F{p)ep 

for all X G hi, where F{p)^ G L^((0,R')), for all g > l,i?' < R, so (fT3H 
becomes 

r F{p)i\p)p<^-Hp=^ 

Jo 

for all radial test functions C,. It is then possible to conclude that F{p)p'^~^ = 0 
a.e. in (0, R), so ([6]) holds and ((p. A) is a solution of (jl]) because of Theorem 
Ol a). 

If ((p. A) is a radial solution of (jl]), say (p(x) = (p{p), D(p{x) = p'{p)ep, then 
(p'(0) = 0. Setting 

hQp{p) 

we have b G C°’"([0, R']) for all 0 < R' < R and 6(0) = 0, so 

fP 2-r' 

■^(p) •= / \b{(j)\^'-^b{a)da 

Jo 

defines a radial function u{x) = u{p) that belongs to C^(f2). An easy com¬ 
putation shows that u satisfies ([8]), hence {u,m,X) is a solution of (jl]) as a 
consequence of Theorem 11.21 bb □ 

Proof of Corollary 11.41 We proceed as in the proof of Corollary 11.31 Let 
(m, m, A) be a solution of ([I]) and set 

i>m'{x) + ho\u'{x)Y ~‘^u\x) m{x) =: F{x) 
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in [a, 6], and F G L‘>{{a,b)). Therefore, 


j-b 

/ F{x)^'{x)dx = 0 

J a 

for all test fnnctions ^ G C°°{[a,b]), so F{x) = 0 for a. e. a: G [a,b]. Hence 
dH]) holds and one implication stated by the corollary follows by Theorem 11.21 
a). Vice-versa, (jS]) holds if we choose u{x) := ~^^b{y)dy, where 

b{y) = {—h'rip'{y))/{hoip{y)) for all y G [a,b], and Theorem 11.21 b) applies. 


□ 
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